Observational Constraints on the Completeness of Space near Astrophysical Objects 
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We consider the observational effects of a deficit angle, w, in the topology of the solar system and 
in the 'double pulsar' system PSR J0737-3039A/B. Using observations of the perihelion precession 
of Mercury, and the gravitational deflection of light due to the Sun, we constrain the magnitude of 
such a deficit angle in the solar system to be 2n{l — w), with < (1 — lu) < 10~® at 95% confidence. 
We calculate the effects of a deficit angle on the periastron advance, geodetic precession rate and 
inclination angle of the double pulsar system and use the observational data to obtain the constraint 
< (1 — ui) < 2.4 X 10~* at 95% confidence. Although this result is weaker than the solar system 
bound, it is in a very different physical environment, where accumulating data is likely to lead to 
tighter constraints in the future. 
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Metric-based gravity theories, like Einstein's and its 
close relatives, are routinely tested using observations 
of astrophysical systems where gravitational fields are 
strong and non-Newtonian. Following Eddington's in- 
troduction of a parameterised form of the metric in 1922 
[l|, and Nordvedt's extension to more general configu- 
rations by including further parameters 0, a sophisti- 
cated frame- work has been devised Q , and a large num- 
ber of careful observations made, in order to place strin- 
gent constraints on the geometry of space-time. However, 
most of these studies have assumed that the astrophysi- 
cal systems they consider should have a trivial space-time 
topology. Under such an assumption, it is possible to 
make strong statements about the magnitude of any de- 
viations from the predictions of general relativity. Here 
we take a different approach: We assume a geometry that 
is locally isometric to the predictions of general relativ- 
ity, but which permits non-trivial global topologies. The 
astrophysical observations can then be used to constrain 
the topology of the systems in question. Of particular 
interest for this study is the recently discovered 'double 
pulsar' system, PSR J0737-3039A/B 

Non-trivial topologies can exist in a number of as- 
trophysically interesting situations, including magnetic 
monopoles, cosmic strings, domain walls and textures @, 
all of which can arise in phase transitions in the early 
universe A tell-tale sign of such topological defects 
is the existence of a 'deficit angle', whereby a wedge of 
space-time appears to have been removed, and the sur- 
faces that remain have been identified. A simple example 
of this feature, first found by Harder in 1959 [st], is the 
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ds^ = -adt^ + + r^de"^ + v?r^ sin^ Odcj)'^, (1) 

a 

where a ~ 1 — 2GM/r and w G (0, 1] is a constant. This 
is well known as the exact solution of apoint-like mass 
on an infinitely thin cosmic string 0, It is clearly 
locally isometric to the Schwarzschild solution, but if we 
allow (j) to run from to 27r, then we see it has a different 
topology. 

The line-element ([I} may initially appear a trivial ma- 
nipulation of the Schwarzschild solution. The effect of 
introducing such a topological defect, however, has non- 
trivial consequences for the geodesic equations. For the 
metric above, these equations have recently been solved 
completely by Hackmann et al. [loj . The missing wedge 
of space-time causes new behaviour that is not present 
in the Schwarzschild solution, such as a precession of the 
ang ular momentum vector about the axis of the string 
By using observational constraints on such phe- 
nomena we can therefore constrain the amplitude of any 
possible deficit angle, and hence place constraints on the 
topologies of astrophysical systems. 

In Section |TT] we summarise how deficit angles affect 
the geodesic equations, and observations of gravitational 
phenomena in the solar system, in Section IHIl we pro- 
ceed to consider the double pulsar. This over-constrained 
system has a large number of well measured observables, 
and provides an excellent laboratory to perform the type 
of test we are considering. In Section IIVI we give our 
conclusions. 



II. SOLAR SYSTEM CONSTRAINTS 

If we choose coordinates so that the orbit we consider 
is in a plane of constant elevation [9 — 7r/2), then the 
effect of a deficit angle on the equations of motion enters 
only through the azimuthal coordinate, 0. To include 
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such a defect, we can then write the range of <{> as 

< (/) < 27r - (2) 

where A(/3 e [0, 2tt) is the deficit angle. Alternatively, we 
can change variables so that /S.ip = 27r(l — w). Now, 
covers the range 

< < 27rw, (3) 
which is equivalent to a simple coordinate redefinition 
(j) w(j). (4) 

A. Time-like orbits 

Using the redefinition ([4]), the geodesic equations be- 
come 



u + w u — 



GM 
IF 



(5) 



where u = l/r, and where primes denote differentiation 
with respect to 0. The angular momentum constant, L, 
is given by 



L^r' sin(6'o) — , 



(6) 



where A parameterises distance along the curve. The 
solution to (HI is then given by 



u = w '^ucRiwc/)), 



(7) 



where ucRi4>) is the usual general relativistic solution 
that occurs when u; = 1. Therefore, we have 

^ = - = [1 + e cos{w(0 - 00 - A0o)}] (8) 

r w L 

where e is the eccentricity of the orbit, 0o specifies the 
phase, and A0o gives the correction due to relativistic 
perihelion precession. 

It can now be seen from ([5]) that the perihelion shift 
per orbital revolution is 



wA0 = 7715— +27r(l-ti;). 



If we now write 



w = 1 — Sw, 



(9) 



(10) 



where Sw G [0, 1) is small, then to first order we should 
expect the precession to be 



A0 ~ A0Gi? + 2tt6w, 



(11) 



where AcjjQn — {A0}^=i is just the usual general rela- 
tivistic prediction with no defect. The precision of the 
agreement with the current data therefore gives a bound 
on Sw^ and hence on w. 



In order to gain observations of A0 from the preces- 
sion of the orbit of Mercury, with respect to the vernal 
equinox of the Sun, it is necessary to take into account 
the precession of the equinoxes on the coordinate sys- 
tem (about 5025" per century), the perturbing effects of 
the other planets (about 531" per century) and the ef- 
fect of the quadrupole moment of the Sun (about 0.025" 
per century), on the perihelion precession of the objects 
that orbit it. These challenges can be addressed in a 
variety of different ways, and the result is authors claim- 
ing slightly different observational bounds on the residual 
A0. Rather than favoring any particular method here, 
we prefer to quote a number of different precise calcula- 
tions, obtained by various authors. These are displayed 
in Table m For further details the reader is referred to the 
original papers, and references therein. For an overview 
of the issues involved the reader is referred to flol. 



Source 



Anderson et al. [14] 
Anderson et al. \15] 
Krasinsky et al. [16] 

EPM1988 

DE200 
Pitjeva [17] 

EPM1988 

DE200 



A(j> - A(j>GR 
(arcsec / century) 



-0.04 ±0.20 
~f0.15±0.14 

-h0.004 ± 0.061 
-0.003 ± 0.061 

-0.017 ±0.052 
-0.011 ±0.052 



10^° X 5w 
(2(7 upper bound) 



6.7 
8.0 

2.3 
2.2 

1.6 
1.7 



TABLE I: The value of the perihelion precession of Mercury 
obtained from observations by various authors, and the re- 
sulting 2(7 upper bound on the permissible values of Sw. The 
acronyms EPM1988 and DE200 refer to different numerical 
ephemerides, which are reviewed in 1 121. We take the sidereal 
period of Mercury to be 0.24 years [l3 |. 

It can be seen from Table U that placing a constraint 
on 5w is not a clear-cut matter, and depends both on 
the data used, and how it is treated. However, all of the 
results in Table U are consistent with the conservative 
statement that the 2a bound on Sw from observations of 
the perihelion precession of IVIercury is within the range 



0<6w< 10" 



B. Null orbits 



(12) 



The geodesic equation for light rays is transformed in a 
similar way, but there is now no L. Under the redefinition 
we therefore have simply 



which has the solution 



(13) 



(14) 
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where, again, ucr{4>) is the usual general relativistic so- 
lution with w = 1. Taking this, together with (f) — tt + S4>, 
we obtain 

1 4GM w5<h TT , , , , 

where R is the impact parameter. The total deflection is 
then 

c, 4GM , , 4GM ^ , , 

wo(p = — h 7r(l — w) = — h TTOw, (16) 



and the deflection caused by the defect is 



PGR 



ttSw, 



(17) 



where Sc/jgr = {50}to=i is the usual general relativistic 
prediction. The bound on Sw, and hence w, can now be 
obtained by comparing ([T7]) with observed deflections. 

The best results available on light bending in gravi- 
tational fields are those of Shapiro, Davis, Lebach and 
Gregory [l^. These authors use almost 2 million obser- 
vations of 541 radio sources by 87 Very-Long-Baseline 
Interferometry (VLBI) sites to calculate the deflection 
caused by the gravitational field of the Sun. Their result 



S(j) = (0.99992 ± 0.00023) S(j)GR, 
which gives us a 2a bound on 6w of 



< (5w < 9.9 X 10 



-10 



(18) 



(19) 



where we have taken the general relativistic prediction 
of light bending for an object whose light ray grazes the 
Sun's limb to be 1.75". This is almost exactly the same 
bound as was achieved for the time-like case I27l. 



III. PULSAR CONSTRAINTS 

The recent discovery of the 'double pulsar' system 
PSR J0737-3039A/B, provides new opportunities 
for testing ideas in relativistic gravity [13, [2l| . These bi- 
nary neutron stars are known to orbit each other with 
a period of 2.45 hours, with much higher velocities 
and accelerations than those found in other binary pul- 
sar systems. By good fortune, this system is also rel- 
atively near to the Sun, and oriented so that we ob- 
serve it nearly edge-on, at an inclination angle of about 
89°. Most importantly, however, and unique to PSR 
J0737-3039A/B, is the feature that both neutron stars 
are detectable as radio pulsars, with periods of 22ms and 
2.7s, respectively for PSR J0737-3039A and PSR J0737- 
3039B. These properties make the double pulsar an excel- 
lent tool for constraining the type of deviant gravitational 
phenomena we are considering in this article. 

To proceed further, we model the space-time geometry 

as 



where 6ij specifies the static geometry of the 3-space. 
We can now assign orders of smallness in the usual way, 
so that the Newtonian potential U is 0(2) small, time 
derivatives add an 0(1) of smallness, and the 3-velocity 
-y* is 0(1) small. 

We will now consider relativistic effects that are good 
candidates for constraining the existence of any deficit 
angle in the double pulsar, before comparing our predic- 
tions with the observational data. 



A. Spin Precession 

Geodetic precession of a body's spin vector about its 
orbital angular momentum vector is a well known predic- 
tion of relativistic gravity, and has already been studied 
in some binary pulsar systems pl| - [26| . However, while 
it took over a decade for geodetic precession to be ob- 



served in the Hulse-Taylor pulsar, PSR B1913-H16 |23|, 

21|. 



it has already been reported in the double pulsar 
We expect this effect to be particularly sensitive to the 
existence of a deficit angle, as it involves integration over 
an orbit. We shall therefore calculate its influence below, 
before proceeding to infer constraints from observational 
data later in the section. 

First, we want the pulsar spin vector, y, to be or- 
thogonal to the world-line of a particle, u^^, and to be 
parallel propagated along that curve, so that 

^^-o^^r n^-r (21) 



where 



dfj.p^'^i^a-, and r is proper time along the 



particle's world-line. The orthogonality condition then 
gives 



Yo = -v% + 0(3), 



(22) 



where v'^ is the 3-velocity of the particle. Multiplying 
through by dr/dt we find that the spatial component of 
Y^^ obeys 

dY 

= To^ov'YJ + Tj.oY^ + Tk^JV^Y'' + 0(5) (23) 

= ~5^jU,oY' + i2~SkjU, + ~6,kU,j - kjU^uyY*' 
+ (1 + 2U)Vk^,v^Y'' + go[fe,.]>^' + 0(5), (24) 



where Tijk = \{^i3,k + 5ik,j 



). The magnitude of 



the 4- vector y should remain constant along the curve 
specified by w^, so that 



|- (g^-'Y^Y,) = 0, 



or, up to 0(2), 



{v%)^ + (1 - 2U)6'^YiYj = constant. 



(25) 



(26) 



ds' 



-(1 - 2U)dt^ + (1 + 2U)Sijdx'dx^ 



(20) 



This implies that a spin 3-vector, 5*', with constant mag- 
nitude, 5^^ SiSj, is given by 

= (1 + U)S, + ]^kjV^v''Sk + 0(4). (27) 
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Alternatively, inverting this expression gives 



= {l-U)Y,-^S,,v^v''Yk + 0{4) 



= (1 + U)S,,Y^ - ^~S^,Skr,yv''Y"' + 0(4), (28) 

which can now be differentiated with respect to and 
dYi/dt can be substituted from above, to give 



dSj 
dt 



where we have used 



dv' 



dt 



For shorthand, one can also write 

D,ki = Sv^j[kU,i] +go[k.i] 

and 

1 



so that 



dS^ 
dt 



0(4). 
(29) 

(30) 

(31) 
(32) 



5'™ + 0(4). (33) 



This is the general expression for spin precession in the 
geometry pO)) . 

To see the effect of a deficit angle on spin precession we 
can now consider a 3-metric that is locally isometric to 
Euclidean 3-space under the coordinate redefinition 
For an elliptic orbit in a plane of constant 0, and with a 
potential of the form U = U{r), we then have — Q and 
Ufi = J7,0 = 0. From ([55t . we then find 



dSr 

~dr 



2 ■ ^ J. ^ 



(34) 



- Ay't'r^y.^u rSr - rw^'t'Sr + -S^. (35) 



dt 



If we now define two new quantities, 5*0 = S^/{r'w) and 
X = {3v'^rU^r/'2' + v'^)w, these equations reduce to 



dSr 

~dt 

dS^ 

"dT 

which have the solutions 



— XS^ 

= — X Sr 



Asini / Xdt 



= Acos<^ / Xdt 



(36) 
(37) 

(38) 
(39) 



where A is a constant of integration. If we now use (|10l) , 
with Sw being 0(2) small, then we find that the preces- 
sion rate integrated over an orbit is 



n = - y Xdt 

3 



27r - ^ / v't'rU^rdt - (f - 5w) / v'^'dt 



GR 



2'k5w, 



(40) 



where VLqh is the general relativistic precession rate, 
when w — Q. The extra precession per orbit due to the 
deficit angle is therefore 2tt5w. 

This result can be seen to be straightforwardly gen- 
eralisable to more complicated orbits if we note that w 
always enters as a multiplicative factor in the 0(2) terms 
in The lowest-order effect of w then occurs when 

it multiplies the 0(0) terms, which are independent of 
gravitational potentials and velocities. 



B. Orbital Inclination 

Binary pulsar observations are sensitive to the sine of 
the inclination of the orbital plane with respect to the 
line of sight of the observer: 



smj 



(41) 



where a is the semi-major axis of the elliptical orbit, and 
X is the 'projected' semi-major axis. Although the value 
of this quantity itself is not dependent on integration over 
an orbit, it is usually determined in terms of another 
quantity that is: The orbital period. 

In the Newtonian 2-body problem the orbital period is 
given by 



(42) 



P,=27T\ —, 

m 



where m = rrip -I- rric is the total mass of the pulsar and 
its companion (another pulsar, in the case of the double 
pulsar) . If there is a deficit angle then we should expect 
this period to be reduced by a factor of w, due to the 
shorter path length, so that 



Pb = 2ttw 



Using this to eliminate a in (j42|) we then find that 
2/3 ( 2'n:\^^^ {nip + m.c)'^^^ 



Pb I ' m 



(43) 



(44) 



where we have also used Xp = xm^/ [nip ~\- m^), which is 
a directly measurable quantity. 
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C. Constraints 

The constraints on periastron advance, inclination an- 
gle and geodetic precession rate obtained in ^20] and [21] 
for the double pulsar system are displayed in Table |lTl 
The first two of these are based on observations made be- 
tween April 2003 and January 2006 using the 64m Parkes 
radio telescope in New South Wales, the 76m Lovell radio 
telescope at the Jodrell Bank Observatory, and the 100m 
Green Bank telescope in West Virginia. The observations 
of geodetic precession were made between December 2003 
and November 2007 using the Green Bank telescope. 



Parameter 


Observed value 


Periastron advance, cj 
Orbital inclination, s = sini 
Geodetic precession rate, Qb 


16.89947 (68) "/year 
0.99974 (-hl6,-39) 
4.77 (-f66,-65) "/year 



TABLE II: Observational constraints on the periastron ad- 
vance, dj, orbital inclination, s — sini, and geodetic precession 
rate of PSR J0737-3039B, Qb, from ,2^] and j21,]. 

Using the expressions found above, together with the 
usual general relativistic {Sw = 0) results, we can now 
write down the periastron advance, inclination angle and 
geodetic precession rates as 



= 3 



5/3 



{nip + rric)'^/'^ 



1--SW 



2n 



2 
2/3 



6w, (45) 



{nip + mc)^/^ 



(46) 



and 



— 



Pj (1 



1 



mc{Anip + 3toc) 



) 2(mp-fm,)4/3 

(47) 

In all of these expressions Pt, xb and e can be mea- 
sured directly, and are given in ^] . The ratio of masses 
iTiA/niB is also known [20j, but not the individual masses 
themselves. We therefore need two observations to ob- 
tain a constraint on the deficit angle. The first will give 



enough information to find the two masses, and the sec- 
ond can be used to compare predicted and observed val- 
ues of a relativistic effect. We choose to use oj as the 
first observable in each case, as it is known to greatest 
accuracy. 

Using the additional constraint from geodetic preces- 
sion, riB, then gives the 2a bound [2^ 



<Sw <1.5x 10" 



(48) 



Alternatively, we can use the inclination angle, s, to get 
the stronger 2a bound 



< Sw <2Ax 10" 



IV. DISCUSSION 



(49) 



We have considered the effects of a deficit angle, w, in 
the space-time metric for the solar system and the dou- 
ble pulsar system. Using observations of the perihelion 
precession of Mercury and the gravitational deflection of 
light we constrain the magnitude of such a deficit angle 
in the solar system to be 27r(l — w), with 



< (1 - w) < 10" 



(50) 



to 95% confidence. Similarly, we have used observations 
of periastron advance and inclination angle in the double 
pulsar system PSR J0737-3039A/B to gain the constraint 



< (1 - w) < 2.4 X 10" 



(51) 



also to 95% confidence. Although this result is weaker 
than the solar system bound, it is in a very different 
physical environment, where accuracy is likely to improve 
in the future, as observational data on the double pulsar 
accumulates. 
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